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Abstract 

We study the cosmological imphcations of the one-loop terms in the string 
expansion. In particular, we find non-singular solutions which interpolate be- 
tween a contracting universe and an expanding universe, and show that these 
solutions provide a mechanism for removing the initial conditions problem pe- 
culiar to spatially closed FRW cosmologies. In addition, we perform numerical 
calculations to show that the non-singular cosmologies do not require a careful 
choice of initial conditions, and estimate the likely magnitude of higher order 
terms in the string expansion. 



^ easther@cfi.waseda.ac.jp 
^ maeda@cfi.waseda.ac.jp 



1 Introduction 



The study of cosmology leads inexorably to an epoch in which the energy den- 
sity of the universe approaches the Planck scale. Superstring theory is the 
leading candidate for a description of physics at the Planck scale, but has not 
developed to the point where the full theory may be employed to construct a 
detailed cosmology. Despite this, the generic properties of superstring models 
can be extracted and their cosmological consequences investigated. In doing 
so, the cosmologist simultaneously hopes to find mechanisms for resolving the 
various problems that remain after the introduction of the inflationary sce- 
nario, and to test string theory in the high energy laboratory provided by the 
early universe. 

Two distinct approaches have been taken to the study of string cosmology. 
Firstly, many cosmologists approximate the full string theory with the first 
terms of the perturbative string expansion IQ, H, H . The non-perturbative fea- 
tures of string theory are implicitly discarded, so this approach breaks down 
at the highest energies. However, since Einstein gravity and conventional par- 
ticle physics agree well with experiment, a successful superstring theory must 
reduce to these theories in the low energy limit. Consequently, it is reasonable 
to expect that at energy scales when stringy effects first become readily ap- 
parent it will be possible to approximate the full string theory by an effective 
theory comprising the usual Einstein- Hilbert action together with higher order 
corrections. In particular, much of the effort to date has focussed upon the 
lowest order, or tree- level action H1-[BT[. 



Gasperini, Veneziano and others have adopted an alternative strategy, 
which they have dubbed the "pre-big-bang scenario" |22, 2^. This is a non- 
singular cosmology, and it attempts to alleviate many of the problems associ- 
ated with the "standard" big bang by appealing to the symmetries of the full 
superstring action, rather than a perturbative limit. The fundamental require- 
ment for such a model is the existence of non-singular or "branch-changing" 
TT| , p!7[| solutions which smoothly interpolate between a contracting universe 



and an expanding one without passing through a singularity. 

The term branch-change originally referred to a non-singular transition be- 
tween two solutions of a simple tree-level string cosmology. The solutions are 



1 



distinguished by a sign choice and are related to one another by a duahty 
transformation. Moreover, this discussion was largely based in the (tree-level) 
string frame, which is related to the usual Einstein frame of general relativity 
by a conformal transformation. A universe which is expanding in the string 
frame may be contracting in the Einstein frame, and it is therefore dangerous 
to simply identify the two branches with "expansion" and "contraction" . How- 
ever, in general the (— ) branch evolves away from a past singularity, while the 
(+) branch represents evolution towards a future singularity [^. The type of 
branch-change required by the pre-big-bang scenario is a non-singular transi- 
tion from the (-I-) to the (— ) branch. The converse, a solution which moves 
from (— ) to (+) is entirely unremarkable, as it corresponds to the generic 
evolution of a spatially closed FRW universe. 

Since string theory is expected to remove, or at least tame, the infinities 
inherent in other models of fundamental physics, it is reasonable to expect 
that it will remove the singularity represented by the big bang in conventional 
cosmology. At tree-level in the superstring action, though, the possibility of 
classical branch- changing solutions to the classical equations of motion has 
been virtually excluded ||2l], although quantum cosmology may pro- 

vide a mechanism for allowing a transition between the two tree-level branches 



26| , p7| . However, the tree- level action is only the lowest order term in the 
full string loop expansion, so consequently in this paper we explore the cosmo- 
logical solutions that arise when one loop contributions from the dilaton and 
modulus terms are included. We show that in this case, the low-energy limit 
of string theory naturally leads to a "bouncing universe", in contrast to the 
tree-level limit. 

Previously, Antoniadis, Rizos and Tamvakis have studied the equa- 
tions of motion derived from the same action with a spatially flat background 
and found some non-singular solutions. However, in this case the scale fac- 
tor increases monotonically from a (non-zero) constant value, which differs 
from the branch- changing solutions envisaged by the pre-big-bang scenario, as 
there is no transition from contraction to expansion. We extend this system 
to the case where the spatial hypersurfaces are allowed to have non-zero cur- 
vature. In doing so we find "bouncing" solutions and investigate the range of 
initial conditions over which they occur. Since we do not directly employ the 
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duality symmetries of the superstring action in our analysis, we do not explic- 
itly identify these bouncing solutions with the branch- changing solutions of 
the pre-big-bang scenario. However, from a phenomenological perspective, a 
bouncing solution resembles a successful branch-change, as in both cases the 
universe is apparently evolving towards a singularity in the distant past and 
away from it in the distant future. 

We also consider the particular initial conditions problem faced by a spa- 
tially closed FRW universe, which in the absence of fine-tuning will typically 



have a lifetime on the order of the Planck scale p9[. This problem persists 
in the presence of inflation: although inflation can begin at the Planck scale 
there is no guarantee that it will do so. Here we show that the contribution 
from the one- loop terms can allow a closed universe to grow arbitrarily large, 
ensuring that the universe will survive long enough for inflation to begin. 

Finally, we numerically integrate the equations of motion for a wide range 
of initial conditions, to demonstrate that the bouncing solutions do not re- 
quire a highly restrictive choice of parameters. We also argue on the basis of 
dimensional analysis that many of the bouncing solutions do not evolve into a 
region where the perturbative expansion is likely to break down, which makes 
it plausible that these solutions will persist if higher loop terms are added to 
the action. 



2 Action and Equations of Motion 



We take as our starting point the one-loop superstring action ||28|, ^ 
1 
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which contains contributions from the Ricci scalar R, the dilaton and a 
modulus field a. The Gauss-Bonnet combination, -R^g, is 
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i?^^"^ - ARf.^R'"' + R^. 



(2) 



Terms which appear at one loop in the string expansion but vanish when the 
metric has the Robertson- Walker form have been dropped from the above 
action. We have adopted the same conventions and notation as Antoniadis, 
Rizos and Tamvakis with the exception of the metric signature, which we have 
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set to (—,+,+,+), and our units correspond to the gravitational constant, 
G = I/Stt. The coefficient A is positive and determined by the four dimensional 
string coupling. The sign of 6 is determined by the relative numbers of chiral, 
vector and spin-3/2 massless supermultiplets, and is proportional to the four 
dimensional trace anomaly of the N = 2 sector of the theory. It will be 
important that 6 can take both positive and negative values. The potential, 
^((t) is defined in terms of the Dedekind r] function, 

e(o-) = In [2e''v\^en] (3) 

where t] is 

oo 
n=l 

Anticipating that ffist and second derivatives of ^ will appear in the equations 
of motion, we note that 

oo 2n7re 
n.=l ^ 

where the subscript denotes differentiation with respect to a. Despite its ap- 
pearance, this is an odd function of a. Furthermore, for large 

^a~-Ysinh(a) (6) 

closely approximates the exact expression, equation (H), and is also antisym- 
metric under a —>■ —a . The accuracy of this approximation for cr ^ could 
be improved by adding terms of the form c„cr" for n = 1,3,5... (since C,a is 
odd, even powers of a will not appear) with the c„ chosen so that the approx- 
imation reproduces the ffist few terms in the Taylor expansion of C,a about 
cr = 0. However, in practice adding these corrections does not alter the qual- 
itative properties of the solutions obtained numerically, so for simplicity we 
have worked exclusively with equation (Q). 

Previously Antoniadis, Rizos and Tamvakis examined the cosmological so- 
lutions for this system that have a spatially flat Robertson Walker metric. We 
extend their work to include the possibility that the spatial hypersurfaces have 
non-zero curvature, and so the appropriate ansatz for the line element is 



1 — kr'^ 



(7) 
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We therefore derive the equations of motion 

3 + A;e-2-) (l + gfoj) - - f = (8) 

2 (cZ. + LU^) (l + 8fu) + (u^ + A:e-2-) (l + 8fu) + + fcr^ = 0, (9) 
4> + 3u^-2§^Rl, = 0, (10) 

^ + 3^^-if^GB = 0, (11) 

where a dot denotes differentiation with respect to t. The Gauss-Bonnet term 
is 

i?L = 24(d) + a;^)(a;^ + Aj (12) 

and / is defined to be 

/ = 1 [Ae^ - Sa<T)] . (13) 

If / vanishes, the system reduces to two free, minimally coupled, scalar fields 
in a Robertson- Walker universe. 

The equations above consist of three second order equations and a con- 
straint, giving a total of five degrees of freedom. In order to facihtate the 
analysis we isolate each of the second derivative terms (remembering that / 
implicitly contains and a) , on the left hand side, giving the following system 

d) = -^^-(cO^ + A.^^, (14) 

= -3cl;0 - 3Ae^ (^cO^ + A j X, (15) 

^ = -3u;& + S^,(u' + ^^ X, (16) 

where 

^ 4 + 2(A0e'^-5(7ea)c^+(cO2 + A;e-2'-)2(3A2e2'^ + 52^2)- ^ ' 

If we wish, we can eliminate any one of the variables by inserting the constraint. 
However, it will be more convenient to work with the equations as they are 
given above, both in the following section when we consider the asymptotic 
form of the solutions and in Section 4, where the constraint will allow us to 
the check the accuracy of our numerical solutions. 
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3 Asymptotic Solutions 



In general this system of equations must be solved numerically. However, con- 
siderable insight into the cosmological properties of this model can be gained 
from analytic considerations alone. In particular, the existence of a bouncing 
solution requires that the scale factor to pass through a (non-zero) minimum 
value and be singularity free. Hence we begin our investigation by considering 
the possible extrema of tu, and a. 

When the scale factor, a = ^ passes through a local minimum, its second 
derivative, a, must be positive. When A; = the constraint, equation (|), with 
cj = requires d = = as well. This is an exact static solution to the 
equations of motion where the values of cj, a and are all arbitrary constants. 
The non-singular solutions of Antoniadis, Rizos and Tamvakis can be regarded 
as the consequence of making a small deviation from this solution in the distant 
past. If the spatial hypersurfaces have negative curvature, k = —\ and the 
constraint cannot be satisfied when ti; = 0. Consequently, the scale factor is 
again monotonic, but the constant solution no longer exists. 

When the spatial hypersurfaces have positive curvature, so that = 1, the 
constraint can be satisfied with cj = and 0, d 7^ 0. Thus, in this case, the 
scale factor may pass through an extremal value. Since the first two terms in 
the denominator of x equal to 4(1 + 8/cj) it follows from the constraint 
that if = 1 the denominator of x is positive, so 

Sign[cZ;|^^o] = Sign[-(8 + ^^'^ - ^d^U)] (18) 

While A is physically restricted to positive values, b can take any real value. 
Since ^aaip^ < for all cr, if 5 > 0, ci) will be negative at an extremum of uj, 
but can take either sign if 5 < 0. This immediately proves the existence of 
solutions for which k = 1 and 6 < where the scale factor possesses a local 
minimum, which is a prerequisite for a successful bounce. However, this does 
not establish the existence of globally non-singular solutions, since it does not 
guarantee that such a solution will always be non-singular. For instance, when 
the tree level action contains contributions from the spatial curvature and an 
axion the Einstein frame scale factor may pass through several local minima, 
but all solutions contain at least one singularity |pT |. 
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Now consider the asymptotic form of solutions to the equations of motion. 
We focus on A; = 1, which is the only choice that can lead to a bouncing 
solution and investigate the properties of a universe which is expanding at 
large positive times 

Type I a(t) > e, t > T, 

Type II a(t) < -e, t> T, (19) 
Type III a{t) ^0, t ^ oo, 

where e and T are both positive numbers. A solution which is contracting at 
large, negative times corresponds to the time-reverse of this case and does not 
need to be examined separately. 



Since the defining condition for inflation is that d > [g^, a solution of 
Type I inflates forever and so the curvature terms will be negligible at late 
times. However, this would imply the existence of solutions for k = that are 
inflationary at late times, in contradiction with the results of Antoniadis, Rizos 
and Tamvakis, so it follows that there are no expanding, asymptotic solutions 
with the form of Type I. Conversely, in the case of Type II, the curvature 
term will dominate the cj^ term where they appear together in the equations 
of motion, and the scale factor will eventually pass through a maximum. Thus 
there are no solutions which expand indefinitely where the scale factor has the 
generic form of Type II. Note that if the scale factor does evolve through a 
local maximum, it may either make a non-singular transition to a subsequent 
stage of expansion or collapse to a singularity. 

Finally, consider possible asymptotic solutions of Type III. Since the Gauss- 
Bonnet term vanishes if d is identically zero, choosing a{t) = a^t is not a 
viable candidate for a late time solution that involves a non-trivial contribution 
from the one-loop terms. Consequently, we need to consider solutions with an 
asymptotic form like a = ^ ait + a2t{hit)'^ which are of Type III but allow 
for the possibility of a non-zero Gauss-Bonnet combination. 

Dropping the contributions from the dilaton terms, make the following 
substitution, where r = In t 

A{t) = ^ = ai + a2r™ (20) 

sit) = ^ = Si + S2T^ (21) 
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which allows us to write equations (|^) and (|TT]) as 

2 




(22) 



(23) 



where A = 7r5/3, a dash denotes differentiation with respect to r and we have 
assumed that cr ^ 1 so that o: . This involves no loss of generality, since 
the full equations of motion are invariant under the transformation a —a. 

With the ansatz, equations ( pOf ) and (^), we can systematically expand in 
powers of r. There are a number of subcases, corresponding to ai,2, si,2, 
and n being either positive or negative, subject to the overall requirement that 
A{t) and S{t) are positive. For k = 1 and 5 < (restricting attention to the 
possible bouncing solutions) we can show that the only asymptotic solution of 
this type is: 

1 



A 



5 

87 



(24) 
(25) 



If <^ \a\ and the scale factor is given approximately by equation (p^) then 
the one-loop terms have a negligible effect in equation (|IOD and it follows that 
in the asymptotic regime the dilaton is effectively constant. As written 
above, this asymptotic form describes an expanding universe as t — oo. A 
universe which is contracting as t — > — oo corresponds to the time reversal of 
this solution. 

We now have the ingredients we need to construct a universe that can 
be arbitrarily large in the distant past, contract to a non-zero minimum size 
and then make a smooth transition to expansion, after which it may grow 
indefinitely. A specific numerical solution of the equations of motion which 
exhibits these properties is shown in Figure 1. The key ingredient in these 
solutions is the presence of the modulus field a and its associated one-loop 
terms. If either 5 > 0, or the one-loop terms are simply absent from the 
action, then the only type of extremum that the scale factor, a(t), can possess 
is a local maximum, and a bounce does not occur. 
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An analogous result can be found for 6 > 0, but in this case no bounce 
is possible. Furthermore, any asymptotic solution to the scale factor-modulus 
system for positive 6 will be equivalent to a similar solution to the scale factor- 
dilaton system in which — oo, although this implies that the system evolves 
into the strong-coupling region and such a solution is therefore unphysical. 

Solutions of this type extend the previous work of Antoniadis, Rizos and 
Tamvakis, further demonstrating that the singularity problem is less acute in 
cosmological models based on the one-loop string action than at tree-level or 
in classical, Einstein gravity. This provides empirical support for the hope that 
a fully non-perturbative string cosmology will resolve the singularity problem 
entirely. 

Physically, the key ingredient in the non-singular solutions is the one-loop 
term that couples the Gauss-Bonnet combination and the modulus; without 
this there are no non-singular solutions. However, the sign of this coupling is 
determined by S and it is only for 5 < that the one-loop modulus terms can 
cause a bounce. Thus, while have estabhshed that adding one-loop terms to 
the low energy string action introduces non-trivial non-singular solutions and 
significantly softens the conclusions reached at tree-level, these corrections do 
not imply a complete absence of singularities. 

The new class of non-singular cosmologies found here shows that string the- 
ory has the ability to resolve the curvature problem typically associated with a 
closed, FRW universe, which has a typical lifetime (and maximum size) on the 
order of the Planck scale. Inflation can easily solve the analogous problem in 
an open universe, which will expand indefinitely, whereas in a spatially closed 
universe inflation is typically only successful if it begins at the Planck scale. 
However, we have shown that the presence of the one-loop modulus terms per- 
mit a closed FRW universe to expand indefinitely. This is not sufficient to 
solve the fiatness problem, but these solutions provide a mechanism for ensur- 
ing that the universe lasts long enough and the energy density becomes low 
enough for infiation to begin without the need for any additional constraints 
in a A; = 1 FRW universe. This is a property of the asymptotic form of the 
solutions, and solutions which contain an initial singularity may still expand 
indefinitely. Furthermore, the analogous asymptotic form for 5 > also al- 
lows indefinite expansion, even though a non-singular, bouncing solution is 
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impossible in this case. 

Solutions of the type displayed in Figure 1 strongly resemble the type of 
branch-change envisaged by the pre-big-bang scenario. Unlike the non-singular 
solutions that exist when k = 0, the scale factor is not monotonic. Ast-^ ±00 
the dilaton is fixed, and the Einstein and string frame scale factors differ only 
by a multiplicative constant. Hence there is no ambiguity in associating a 
solution which is contracting in the string frame at early times to one that is 
contracting in the Einstein frame, as there might be if the dilaton were evolving 
rapidly. 

We can also gain some insight into the behavior of the dilaton and modulus 
fields. Consider equations (p!5| ) and (|I6|), and recall that the denominator of x 
is positive for /c = 1, so 

Sign[0|^J = Sign[-{8-6&%^)], (26) 
Sign[a|_o] = Sign [^^.(a)] , (27) 

since A > on physical grounds. Thus when 6 < 0, which is the case we are 
most interested in, the dilaton can possess both maxima and minima and, at 
least in principle, may undergo several oscillations before the solution becomes 
established in one of the asymptotic regimes. However, any extremal value of 
a when a < is a local maximum, and any extremal value when a > is a 
local minimum. Thus when 6 < 0, a can have at most one extremum as it 
evolves from one asymptotic regime to the other. 

When the scale factor a passes through its minimum value, a > 0. As this 
is the minimal requirement for the existence of inflation, it follows that in a 
pedantic sense a bouncing solution is also inflationary. Obviously, the astro- 
physical constraints that a successful inflationary model must satisfy are much 
more demanding than simply requiring that the scale factor undergo positive 
acceleration. At late times the asymptotic form of the solution approaches 
the borderline condition between inflationary and regular growth, as a — > 0, 
which resembles the coasting solutions that have been discussed in the context 
of standard inflationary cosmology pi . 



In the next section we use numerical techniques to investigate the generality 
of these non-singular solutions. However, note that while the scale factor is 
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monotonic if = — 1, there is a simple exact solution 

uj = cjo + lnt, (28) 
= 00, (29) 
a = (To, (30) 

which describes an empty, curvature dominated universe, and a = 0. However, 
unlike the linearly expanding solution found for k = 1, this solution holds when 
A and 6 are zero and the Gauss-Bonnet combination will vanish exactly, so it 
is not related to the one-loop terms in the string action. 



4 Numerical Results 

In the previous section we displayed a particular example of a non-singular 
universe, associated with the one-loop terms in the perturbative expansion 
of the superstring action. In this section we estimate of the likely impact of 
higher loop terms and investigate the range of initial conditions that give rise 
to a non-singular universe. 

While it is simple to determine which solutions are singular and which are 
non-singular, the more qualitative distinction between a solution which comes 
"close" to the Planck scale, and one where the bounce occurs at a considerably 
lower energy scale is obviously important. One way to quantify this is through 
the maximum values of the scalar quantities that can be constructed from 
the metric curvature, such as R, R^yR^^ , R^^kxR'^'^'^^ or various higher order 



combinations. In particular, Brandenberger, Mukhanov and others |3^, 
have constructed models that ensure these quantities remain sub-Planckian for 
all homogenous and isotropic cosmological solutions. The action considered by 
us does not have this property since it does admit some singular solutions, for 
which the curvature invariants will exceed any given finite value. However, 
the curvature invariants derived from the non-singular solutions are typically 
sub-Planckian. As a specific example, the value of R^^i^xR^^'^^ for the solution 
of Figure 1 is plotted in Figure 2. 

While it is reassuring that the curvature terms do not exceed the Planck 
scale, the dilaton and modulus fields couple directly to higher-order curvature 
terms, and if the impact of higher terms in the perturbative expansion is to 
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be small then these combinations must be sub-Planckian as well. Without a 
detailed calculation at two-loop order and beyond, it is impossible to state 
definitively whether the non-singular solutions found at one-loop will persist if 
higher order terms are added to the action. However, on dimensional grounds 
the two-loop terms can be expected to be smaller than the one-loop Gauss- 
Bonnet term provided that u, and a do not become significantly greater 
than unity (in Planckian units) while a = e"^ remains less than unity. Also, if 
e"^ (which is effectively the loop expansion parameter) becomes large, we can 
expect a significant contribution from higher-order terms. Hence on the basis 
of dimensional analysis alone we expect that the two-loop terms will be at 
worst roughly equal to the one-loop terms if 

dj, a, e't', 1/a < 1 (31) 

at all times. Note that this is typically a harsher constraint than requiring 
that the curvature invariants do not exceed the Planck scale. 

Having established a heuristic criterion for the importance of higher order 
terms, we also wish to establish whether bouncing solutions require a special 
choice of initial conditions. Given A and 6, a solution to the equations of 
motion is fully specified by five initial conditions, the sixth being fixed by 
the constraint. The set of points in "initial conditions space" for which the 
system evolves towards a particular attractor (that is, a particular asymptotic 
solution) is known as the "basin of attraction" of the corresponding attractor. 
The full system of equations is five dimensional, and therefore cannot be easily 
visualized using a phase space approach. Since the non-singular dynamics 
depends on the behavior of the modulus field we could set A to zero (or drop 
the dilaton entirely) and reduce the system to three dimensions. Previously, 
however, Cornish and Levin |^ have examined the basins of attraction for 
two field infiationary models, whose equations of motion are similar to those 
considered here, by numerically integrating to find the asymptotic behavior 
for many choices of initial conditions. Using this technique we do not have to 
restrict ourselves to a reduced system in which the dilaton is trivial, and we 
can investigate the basin of attraction for solutions where u) and a have the 
asymptotic form of equations ( [2D| ) and ( PT| ) and the dilaton tends towards a 
constant. 
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In Figure 3 we display the asymptotic form of the solutions derived from a 
sequence of two dimensional slices through the initial conditions space, where 
each slice represents a 400 x 400 grid of c^o and (Tq values, for five choices 
of (Jo between and 2. Since the equations of motion are unchanged by the 
transformation a — cr there is no need to carry out separate integrations for 
(Jo < 0. The integrations were carried out using the Bulirsch-Stoer method 
^ , with the form of the equations of motion given by equations (0) to ([TBI) , 
while the constraint was used to check the accuracy of the numerical routines. 

The intersection of the basins of attraction for a universe that is non- 
singular as t ^ ±oo is shown in Figure 3. The intersection of these two 
attractors defines the region of initial conditions space for which leads to a 
non-singular, bouncing cosmology. Since this is a substantial volume of the 
total range of initial conditions, it follows that the universe does not have to 
be "fine tuned" in order to ensure that it is non-singular. Furthermore, the 
results obtained here do not depend strongly on 0o or 0o, or the magnitudes 
A and 5. 

As well as displaying the choices of initial data corresponding to a bouncing 
solution. Figure 3 shows the subset of points for which the higher-loop terms 
in the perturbative expansion are not expected to make a substantial contri- 
bution. We see that this in turn represents a large subset of the total basin 
of attraction, suggesting that the existence of non-singular solutions is not a 
quirk of the one-loop action. Thus it is clear that the non-singular solutions 
are in no way special, and that a substantial subset of them are such that the 
values of a, 0, a and their corresponding velocities are sufficiently small at 
all times to ensure that the action is unlikely to be dominated by higher-loop 
corrections. 



5 Discussion 

In this paper we present explicit non-singular cosmological solutions derived 
from the one-loop superstring action. These solutions represent a FRW space- 
time where the scale factor makes a transition from contraction to expansion, 
while remaining non-zero at all times. As such, they resemble the type of 
branch-changing solution that is a prerequisite for the successful implemen- 
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tation of the pre-big-bang scenario, and demonstrate that ahhough branch- 
changing does not occur at tree level in the superstring action, non-singular 
"bounce" solutions do exist at one-loop. 

While the non-singular solutions described here occur only when the spa- 
tial hypersurfaces have positive curvature and the parameter 6 in the one-loop 
action is negative, our numerical calculations establish that the non-singular 
solutions do not require a careful choice of initial values for the scale factor, 
dilaton and modulus fields. Furthermore, while we have not made a detailed 
analysis of the two-loop terms we expect that their contribution will not dom- 
inate the one-loop terms, providing cause for cautious optimism that these 
non-singular solutions will not vanish when higher loop terms are incorporated 
into the action. 

These solutions allow us to address the initial conditions problem peculiar 
to the k — 1 FRW universe, which typically has a lifetime on the order of the 
Planck scale. While GUT-scale inflation will solve the other difficulties faced 
by the standard model of the big bang, inflation can only prevent a closed FRW 
universe from recoUapsing if it commences at the Planck scale. However, we 
have seen here that for a wide range of initial conditions the equations of 
motion derived from the one-loop string action permit a closed FRW universe 
to expand indefinitely. Consequently, this provides a string motivated solution 
to the version of the curvature problem faced by a A; = 1 universe without 
requiring that inflation commences at the Planck scale. 

The inflationary epoch cannot be described within the context of the model 
discussed here, as we have not included matter fields in the action that can 
drive inflation followed by a graceful exit to a non-inflationary universe with 
the matter content we observe at low energies. Thus, while our solutions 
describe a universe which grows arbitrarily large, we know that in practice 
the assumptions that underpin the action we consider must break down at 
sufficiently low energy scales. Furthermore, the existence of solutions that 
describe a closed FRW universe which grows arbitrarily large require that the 
modulus field is also evolving continuously, which is made possible by the 
fact the modulus field contains no explicit mass terms at tree level. In a 
realistic theory the moduli fields must become massive after supersymmetry 
breaking, which puts a lower limit on the validity of our action even in the 
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absence of other matter fields and non-perturbative effects. Thus attempts 
to use the properties of these solutions in a more realistic cosmological model 
must provide a mechanism that will allow a "graceful exit" from the modulus 
dominated expansion into an inflationary phase. Furthermore, the late time 
evolution of the dilaton must be small enough to ensure that the constraints on 
the time variation of G, the Newtonian gravitational constant, are obeyed. We 
do not address this problem here, but note that this difficulty afflicts almost all 
string inspired models, and is not peculiar to the specific example considered 
here. 

Physically, we have seen that the key difference between the one-loop and 
tree-level actions that permits the existence of non-singular solutions is the cou- 
pling between the modulus and the Gauss-Bonnet combination. However, the 
presence of this term alone cannot guarantee a non-singular universe. Firstly, 
while the parameter S can take on both positive and negative values, non- 
singular solutions only arise when 6 < 0. Secondly, even if 6 is negative, non- 
singular evolution is not guaranteed. However, the existence of non-singular 
solutions to the one-loop action represents a significant improvement on the 
results found at tree-level and shows that the perturbative limit of string the- 
ory can address the singularity problem that is characteristic of the standard 
big-bang. 

Generically, extensions to Einstein gravity consisting of combinations of 
scalar fields and higher order curvature terms will not lead to a singularity- 
free theory, but are more likely to make the existing singularity problems worse. 
Thus it is noteworthy that string theory naturally leads to a low-energy action 
in which non-singular solutions are possible. More general theories have been 
considered by Brandenberger and Mukhanov et al. who derive a higher 

order gravitational action for which all homogeneous and isotropic solutions 
are non-singular. In doing so, they show that this requirement places a strong 
constraint on the possible form of the action, in the absence of fundamental 
scalar fields. If fundamental scalar fields were incorporated into their approach, 
it would be possible to investigate whether models containing second order 
curvature invariants (such as the Gauss-Bonnet combination) coupled to scalar 
fields generically permit bouncing solutions and to what extent this is a special 
property of the one-loop superstring action. Conversely, Rizos and Tamvakis 
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p9| show that with for a scalar field coupled to the Gauss-Bonnet term with 



a spatially fiat metric, non-singular solutions are possible for a wider range 
of couplings than the function ^(a) that is derived from string theory, and it 
would be of interest to extend their analysis to the k = 1 case. 

The fact that incorporating string loop contributions into the action leads 
to a non-singular cosmology is reminiscent of previous attempts to solve the 
singularity problem by including higher-order terms derived from quantum 
corrections to the gravitational action PO, HT], E2[ , which have long been known 



to remove, or at least soften, the singularity associated with the standard model 
of the big bang. 

Superstring theory currently offers us a tantalizing glimpse of a paradigm 
that will provide a unified description of fundamental physics which modifies 
the predictions of general relativity in such a way as to remove the initial 
singularity associated with the big bang. At present it is impossible to know 
whether this promise will be fulfilled, but there are strong indications that 
string theory addresses most, if not all, of the problems of the conventional 
big bang that can be traced back to the Planck scale. 
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Figure Captions 



Figure 1: A particular non-singular solution is displayed, where 6 = —16 x S/vr 
and A = 1. The chosen initial values are (Tq = 0, 0o = ^3, 0o = —5 x 
10~^, 6"o = 0.08 and = 0.01, so that the constraint requires Uq = 3.25114, 
and all quantities are expressed in Planckian units. The scale factor passes 
through a minimum value of approximately 10, and the solution approaches 
the asymptotic form of equations (p^ and p5| ) as t — > ±oo. 



Figure 2: The value of the scalar, R^y^\R^" is plotted for the solution de- 
picted in Fig. 1, demonstrating that it remains well below the Planck scale at 
all times. Note the different scale on the time axis in this plot. 



Figure 3: Five different slices through the initial conditions space are displayed, 
with 00 = —3, 00 = —0.005, A = 1 and 5 = —16 x 3/tt. Note that not all 
choices of initial conditions correspond to a universe with k = 1. The set 
of initial conditions leading to a non-singular universe constitutes the basin 
of attraction for the bouncing solutions. The regions inside the solid lines 
correspond to solutions which do not violate the approximate bounds on the 
validity of the one- loop approximation, equation (|3lD . 
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